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I. INTRODUCTION 


Natural gas as a source of clean and high-quality 
energy has enormous potential for use as a type of energy 
resource. With the development of the natural gas industry 
and the establishment of the gas pipelines and distribution 
networks, situations of possible leaks have been shown to 
be a major problem in economic matters and, mainly, in 
safety. Any leaks in natural gas pipelines will certainly 
cause contamination and environmental pollution 
problems. Threats become even more serious when natural 
gas and air mix in explosive concentrations, which may 
cause fires and even explosions. 


The study of leaks in distribution networks and natural 
gas pipelines depends directly on the rate of leakage that 
has occurred, which will be the basis for studying the 
diffusion of gas and for predicting the consequences of the 
accident. In this line, it is possible to highlight the work of 
Montiel (1998), in which he proposed the study of models 
of holes with large dimensions, from holes of small 
dimensions. Zou et al. (2008) analyzed the thermal process 
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Abstract — The natural gas leakage rate plays an important factor in 
simulating leaks in gas transportation and distribution pipelines. The 
method for leakage rate calculation involves identifying the geometry and 
type of leak opening, associated with the application of the theories of 
thermodynamics and gas dynamics, as well as energy conservation 
equations. In this article, a review and analysis of the flow of gas released 
into the atmosphere in cases of accidents of this nature is carried out. 
Such a situation is critical, given that certain leaks can create very 
extensive danger zones around the immediate vicinity of the leak point. 


of emission from a natural gas reservoir, establishing a 
corresponding mathematical model. Azman et al. (2011) a 
leak model established for the completely ruptured pipes, 
analyzing the corresponding pressure drop along the 
pipeline. 

Hou et al. (2011) proposed the use of the FLUENT 
software, which corresponds to a computer simulation 
software in the area of fluid dynamics, aiming to simulate 
the diffusion of natural gas in the atmosphere as a result of 
leaks. 


The process of leakage and diffusion of natural gas into 
the atmosphere characterizes the condition of a diffuse 
accident, with no reliable basis for an effective solution. In 
this context, the initial forecast of the volume of gas 
released into the atmosphere, after the occurrence of the 
leak, would allow initial actions to quantify the incident 
and, consequently, take appropriate actions to reduce risks 
and losses. 
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Il. THE NATURAL GAS LEAKAGE RATE IN 
PIPELINES AND DISTRIBUTION NETWORKS 


European Pipeline Incident Data Group (EGIG) 

classifies pipeline leaks into three distinct groups: a) leaks 
occurring from small holes (leak hole diameter less than 
20 mm); b) leaks occurring from holes with large 
dimensions (hole diameter greater than 20 mm, but smaller 
than the pipe diameter); c) complete rupture of the 
pipeline. 
2.1 — Leaks occurring from small holes - for the present 
case, the following simplifying assumptions are assumed: 
(a) the flow is isentropic at the gas release point and 
adiabatic in the pipeline; (b) the flow model is essentially 
1D; (c) the pipeline can be regarded as non-deformable 
and the interaction between the fluid and the pipeline due 
to vibration problems is negligible. From the differential 
equation of Euler kinematics, the adiabatic equation and 
the equation of an ideal gas, we have: 


Fig. 1: Schematic of possible leakage 


2k ZRT, 
U3 = Ta 


k-1M 


E 


valid, from the condition: u2 << us. 


Replacing Eq. (1) in the leakage mass flow equation, we 
can get the leakage rate: 


Dcp al 2K M 
Q = Cop2 A |E IZRT, 


j (2) 


e-o” 


The classification of the flow at the leakage point as 
sonic or subsonic can also be performed by comparing it 
with the so-called critical pressure ratio (CPR) given by: 

k 
CPR ae m = 
7 (; + 


given: Q = mass flow [kg/s]; 


A = hole area [m?]; 
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M = molecular weight [kg.kmol*']; 
p = gas pressure [Pa]; 

Patm = atmospheric pressure [Pa]; 
T = gas temperature [K]; 

Z = gas compressibility factor; 

u = gas velocity [m/s]; 

Co = discharge coefficient (= 0.72); 
k = specific heat ratio; 


R = 8314 = universal gas constant [J.kmol!.K"!]. 


The ratio (Patm/p2 > CPR) corresponds to a subsonic 
flow while the ratio (pam/p2 < CPR) corresponds to a 
sonic flow. For the sonic flow condition, equation (3) can 
be written in the following simplified form: 


k+1 (4) 
EP T Mk ( 2 ‘a 
Q = Copa A loan \k +1 
and, for subsonic flow: 
Q = Cop2A 2k- M (mý - (Pern) ® ©) 
= topat Leal ZED Dy 
2.2 — Leaks occurring from large holes - in these 


circumstances it becomes necessary to evaluate the flow 
conditions existing in the hole and in the pipeline. 


2.2.1 - Subsonic flow in the tube and sonic flow in the hole 
- if confirmed, the following conditions: 


(2) < CPR 
p2 


6 
SAN 2 + (k — 1)Ma,? (6) 
P2 > Pima, k41 


Eq.(4) prevails for calculating the rate of gas released into 
the atmosphere. 


Ma corresponds to the Mach number, whose value 
translates a ratio between the relative speed of the gas and 
the speed of sound in this medium. Your expression is: 


Ma = £ (7) 
KRT 


2.2.2 - Subsonic flow in the pipe and at the hole — now if 
confirmed, the following conditions: 


(2) > CPR 


P2 


8 
Sah 2 + (k — 1)Ma,’ (8) 
P2 7 Pima, as 
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Eq.(5) prevails for calculating the rate of gas released into 


the atmosphere. 


2.3 — Leaks that occur in the complete rupture of the 
pipeline - when the pipeline is completely ruptured, the 
mass flow rate of the natural gas pipeline: 


Q = Cop2A 
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2M k 
ZRk-1 Ey E E 


p2 


T-T 


Pı 


(9) 


II. PROPERTIES OF NATURAL GAS 


Considering a typical chemical composition for natural 
gas - natural gas distributed by PETROBRAS in southern 


Brazil - we have: 


Table.1: Natural Gas — Typical composition chemical 


Components Molar Molar Pc T: (K) 

fraction mass (MPa) 

(%) (kg/kmol) 

CH4 0.8901 16.043 | 4.596 | 190.6 
C2H6 0.0593 30.069 | 4.883 | 305.4 
C3Hs 0.0185 44.096 | 4.250 | 369.8 
n-C4Hio 0.0042 58.123 3.796 | 425.2 
i-C4H 0 0.0031 58.123 3.648 | 408.2 
n-CsHi2 0.0011 72.151 3.370 | 469.70 
i-C5H12 0.0008 72.151 3.380 | 460.4 
n-C6H14 0.0008 86.178 3.010 | 507.5 
N2 0.0067 28.013 3.400 | 126.3 
CO2 0.0154 44.010 | 7.382 | 304.2 
Natural Gas | 1.0000 18.374 | 4.632 | 204.48 


Table.2: Natural Gas — Specific heats [J/kg.K] 


Components Molar Cp Cy 
fraction 
(%) 
CHa 0.8901 | 2253.3 | 1735.1 
C2H6 0.0593 | 1754.3 | 1477.8 
C3Hs 0.0185 | 1672.9 | 1484.4 
n-C4H10 0.0042 | 1708.9 | 1565.9 
i-C4H10 0.0031 | 1673.5 | 1530.5 
n-C5H12 0.0011 | 1649.1 | 1533.9 
i-C5H12 0.0008 | 1645.0 | 1529.8 
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n-C6H14 0.0008 | 1673.1 | 1576.6 
N2 0.0067 | 1037.5 | 740.71 

CO2 0.0154 | 869.34 | 680.43 
Natural Gas | 1.0000 | 2096.0 | 1643.6 


C, 2096 (10) 
i Cy 1643.6 | eee 
Z = 1,00 (11) 


IV. LEAKS - CASE STUDY 
4.1 — Leaks occurring from small holes — Case study 


A section of carbon steel pipe with a diameter of 8 in 
(= 203.2 mm) is also considered, transporting natural gas 
at a pressure pı = 35 bar = 3.5 MPa, being subjected to a 
leak located at a point (L) where the pressure p2 = 27 bar = 
2.7 MPa. For a supposedly constant gas flow temperature 
of 25°C, we have (Pam = 1 bar = 0.1 MPa): 


k 1.275 


2 \k-1 2 0.275 
(ei a 


k+1 


(Pe) < CPR — sonic flow 


P2 


Q = Cop2 A 


2.275 


(18.374)1.275 


2 \o275 
Q = (0.72)(2.7) A aienq ar) 


Q = (3.50(10-%))A = (2.75(1073)) a? 


Assuming holes varying between 2 and 20 mm (small 
holes), the values are obtained: 
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Table.3: Hole diameters x flow of gas released into the 


atmosphere 
Hole diameters (mm) | Leakage mass flow (kg/s) 
2 1.10 x 10° 
4 4.40 x 10% 
6 9.90 x 108 
8 1.76 x 107 
10 2.75 x 107 
12 3.96 x 107 
14 5.39 x 107 
16 7.04 x 107 
18 8.91 x 107 
20 1.10 x 10° 


Hole diameters (mm) x flow of gas (kg/s) released 
into the atmosphere 


2 4 6 8 10 12 14 16 18 


Fig. 2: Hole diameters x flow of gas released into the 
atmosphere 


4.2 — Leaks occurring from large holes — Case study 


Keeping the same flow assumptions defined in the 
previous item and a velocity u; = 20 m/s, we have: 


k 1.275 
ron 2 k-1 2 0.275 0.55 
7 G + 7) (zz) E 
0.1 
Patm _ — = 0.038 
P2 2.7 
(=) < CPR 
P2 
Uy 20 
Ma, = = —__________ = 0,048 


KRT, 1.275(8314)298.15 
M 18.374 
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2 + (k — 1)Ma;? 


>p,M 
P2 7 Pima, k41 


2 + (0.275)0.048? 


p2 > 3.5(0.048) — OK! 


2.275 
Thus: 
k+1 
ae r a 2 j~ 
Q = Copa A IZRT,\k F1 


2.275 


(18.374)1.275 ( 2 j” 


EE A (1)8314(298.15) \2.275 


Q = (3.50 (107?))A = (2.75(107?))d? 


Assuming now holes varying between 20 and 200 mm 
(smaller than the pipe diameter), the values are obtained: 


Table.4: Hole diameters x flow of gas released into the 


atmosphere 
Hole diameters (mm) | Leakage mass flow (kg/s) 
20 1.10 x 10° 
40 4.40 x 10° 
60 9.90 x 10° 
80 1.76 x 10° 
100 2.75 x 10° 
120 3.96 x 10° 
140 5.39 x 10° 
160 7.04 x 10° 
180 8.91 x 10° 
200 1.10 x 104 
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Hole diameters (mm) x flow of gas (kg/s) released 
into the atmosphere 


Fig. 3: Hole diameters x flow of gas released into the 
atmosphere 


V. CONCLUDING REMARKS 


This article presents a comparative study of the volume 
of natural gas released into the atmosphere in situations of 
eventual leaks in the pipeline. In addition to the 
mathematical review of the respective models, a 
comparative analysis of problems is carried out 
considering small holes (< 20 mm) and holes with larger 
dimensions (above 20 mm, but below the pipe diameter). It 
is noted that the flow of gas released into the atmosphere is 
not negligible, identifying that the problem must be dealt 
with responsibly and in a time of action as soon as 
possible. Obviously, in these circumstances, it is assumed 
that the eventual leak is easily detectable by one of the 
several alternatives available for monitoring this type of 
situation. 
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